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摘 要 : 本 文 使 用 快速 鲁 棒 性 主 成 分 分 析 (Fast RPCA, Fast Robust Principal Component 
Analysis) 方 法 对 日 昂 序 列 图 像 中 的 日 网 喷 流 活动 进行 检测 。 检 测 的 基本 思路 是 利用 Fast 
RPCA ee 硫 分 解 的 思想 与 日 网 序列 图 像 中 有 着 变化 尺度 稍 小 且 占 比较 大 的 随 
机 变化 背 最 成 分 、 变 化 尺度 较 大 且 占 比较 小 的 日 网 喷 流 的 特点 相 结 合 ， 以 此 来 实现 随机 复杂 
ZEW ae RoR SRE SD 目标 之 间 的 分 离 ， Aire bE A ATR RCW A RW. KAK 
阳 动 力学 天 文 台 (Solar Dynamics Observatory, SDO) I £ Ł K ^ a 4% 28 44 (Atmospheric 
Imaging Assembly, AIA) 观测 设备 的 两 组 不 同时 间 段 、 不 同 波段 、 不 同 观测 位 置 的 日 网 序 
列 图 像 作为 研究 对 象 。 研究 内 容 主要 包括 日 网 序列 图 像 的 预 处 理 、 gl Fit ADM Fast RPCA 

1 与 帧 间 差 分 法 的 检测 结果 对 比分 析 。 实 验 结果 表明 ， 与 帧 间 差 分 法 相 比 ，Fast RPCA 7 

能 够 检测 出 强度 较 弱 的 日 晃 喷 流 且 提 高 了 日 园 喷 流 检 测 的 准确 度 。 
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磁力 线 运 动 的 等 离子 体 流 , 在 不 同位 置 和 不 同 波段 观测 到 的 喷 流 具有 不 同 的 形态 学 和 动态 学 
特征 后 。 近 年 来 的 研究 表明 ， 日 铣 喷 流 在 分 析 太 阳 磁 场 拓扑 结构 ， 加 热 局 地 日 晃 以 及 加 速 太 
阳 风 等 方面 可 能 起 着 十 分 重要 的 作用 ,同时 喷 流 活动 经 常会 伴随 着 耀 斑 、 暗 条 和 日 冕 物质 抛 
射 等 活动 现象 Bl 

文中 通过 随机 选择 20 ^X BY BBCI A Bs IB: Extreme Ultraviolet, EUV) Siji a Ht 
可 知 ， 这 20 个 喷 流 都 起 源 于 迷你 浊 条 爆发 ， 并 且 上 暗 条 的 喷发 始 于 喷 流 基部 亮点 位 置 附近 。 
在 喷 流 喷发 过 程 中 存在 迷你 暗 条 的 上 升 运动 ,在 暗 条 到 达 喷 流 顶 端 位 置 之 前 , 测量 得 到 的 暗 
ei X 射线 喷 流 的 平均 宽度 ， 这 与 喷 流 喷发 是 由 微 暗 条 喷发 驱动 的 观点 一 致 。 
文 辐 首次 对 日 铣 爆 裂 喷 流 ， 和 由 它 引起 的 同时 发 生气 泡 状 和 喷 流 状 CME 的 活动 现象 进行 分 
i. He tH SAN ARE H Sa PALING AA BY Oe PS ETAL PA AP BA Oa YNZ AS HK A 
CME 现象 进行 解释 : 喷 流 冷 分 量 是 由 喷 流 底部 的 暗 条 喷发 引起 , 并 进一步 导致 气泡 状 CME; 
喷 流 状 CME 与 喷 流 热 分 量 有 关 ， 热 分 量 是 限制 暗 条 的 闭合 场 与 周围 开放 场 磁 力 线 重 新 联接 
而 产生 的 向 外 移动 的 热 等 离子 体 。 文 中 使 用 SDO 卫星 上 的 AIA、 日 震 学 与 磁场 成 像 仪 
(Helioseismic and Magnetic Inager，HMD 观 测 设 备 的 图 像 数 据 和 lm 新 真空 太阳 望远镜 (New 
Vacuum Solar Telescope，NVST) 的 Ho 线 心 图 像 数据 对 喷 流 的 演化 和 喷 流 源 区 两 个 相 邻 暗 条 
丝 间 的 变化 进行 分 析 发 现 , 该 双向 环 状 喷 流 的 爆发 是 由 于 喷 流 源 区 相 邻 暗 条 细 丝 间 的 磁 习 
联 所 导致 ， 和 目前 普遍 认为 的 双向 喷 流 的 触发 机 制 有 很 大 的 不 同 。 
随 着 太阳 观测 设备 的 发 展 , 高 时 空 分 辩 紊 的 日 罗 图 像 数 据 越 来 越 多 , 如 何 从 海量 的 数据 
中 ， 通 过 图 像 处 理 技术 ， 将 日 冕 喷 流 自动 检测 出 来 ， 将 有 助 于 研究 人 员 对 喷 流 的 形成 过 程 、 
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FAC AL til] DA De PBT RY H S8 KADIN MAN A BE UJI Ss E Ja ES EE PE OE EP SY BB A) TE 
流 与 暗 条 、CME 等 伴随 活动 现象 之 间 的 内 在 联系 提供 便利 。 
se a Sobel 边缘 检测 方法 和 Hough 变换 相 结合 来 对 SDO/AIA WLW ACHE HAY H et 
进行 检测 ,并 将 检测 出 的 日 晃 喷 流 相 应 的 活动 持续 时 间 , 长 宽度 和 视线 速度 等 参数 存 入 日 球 
物理 学 事件 知识 库 (Heliophysics Event Knowledgebase，HEK) 中 。Hough 变换 将 检测 对 象 从 
目标 空间 转换 到 参数 空间 ,避免 了 在 目标 空间 检测 时 的 目标 分 类 运算 ,但 是 该 方法 比较 依赖 原 
始 图 像 空间 中 检测 对 象 的 曲线 表达 式 。 而 日 园 喷 流 在 形态 ;上 是 复杂 多 变 的 ， 其 边缘 曲线 的 表 
达 式 会 难以 确定 ， 造 成 检测 结果 不 佳 。 文 &9] 中 使 用 帧 间 差 分 法 对 日 出 (Hinode) 卫 星 上 X HY 
线 望 远 镜 (X-ray Telescope， a ara E 并 对 检测 到 的 喷 流 的 长 度 、 速 度 等 
参数 进行 计算 和 统计 分 析 。 帧 间 差 分 法 原理 简单 、 易 于 实现 ,但 是 在 面 对 喷 流 运动 速度 缓慢 
和 强度 较 弱 时 检测 结果 不 理想 。 文 外 中 首先 对 NVST 观测 到 的 多 个 Ha 日 浪 事 件 进行 观测 研 
究 ， 得 到 该 类 事件 的 观测 特性 和 主观 判断 结果 ， 然 后 对 观测 数据 上 比较 暗 的 结构 (黑子 、 上 暗 
条 和 和 针 状 体 等 ) 进行 标记 ,在 利用 阔 值 分 割 的 方法 分 离 出 观测 数据 中 的 上 暗 结 构 ， 排 除 已 经 标 
Pa 将 剩余 的 暗 结构 和 主观 判断 结果 进行 比较 ， 最 终 得 到 Ho 日 浪 。 该 方法 需要 对 
黑子 ， 暗 条 等 结构 标记 ， 比 较 费 时 。 

FE, 日冕 喷 流 的 检测 可 以 看 成 是 在 日 冕 序列 图 像 中 检测 突变 , 或 者 是 说 从 复杂 随机 
7 变化 的 背景 中 将 有 规律 变化 的 喷 流 活动 检测 出 来 。 因此, WAT SE ELS RR eS it BT E 
> 的 分 离 是 检测 日 冕 喷 流 的 关键 。 鲁 棒 性 主 成 分 分 析 方 法 在 运动 目标 检测 中 有 着 广泛 的 应 用 。 
本 文 用 该 方法 对 SDO/AIA 观测 数据 中 不 同时 间 段 、 不 同 波段 、 不 同 观 测 位置 的 日 冕 喷 流 进 
行 检 测 。 检测 结 果 表 明 该 方法 与 传统 的 帧 间 差 分 法 相 比 能 够 检测 出 更 加 完整 的 喷 流 ,检测 出 
的 喷 流 区 域 孔洞 较 少 、 比 较 连续 。 

本 文 将 从 以 下 四 个 部 分 对 所 做 工作 进行 展开 和 总 结 : 1. 对 所 使 用 的 两 组 数据 的 特点 和 差 
异 进 行 描述 ;2.Fast RPCA 方法 的 原理 及 其 应 用 于 喷 流 检测 的 可 行 性 分 析 ; 3. 基于 Fast RPCA 
的 日 晃 喷 流 自 动 检测 结果 及 相关 参数 分 析 。4. 结 论 与 存在 的 问题 。 
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1 SDO/AIA 实验 数据 


SDO 是 美国 太空 总 署 (National Aeronautics and Space Administration, NASA) F 2010 年 2 
J Beat - 28 A, KER AIA 观测 设备 成 像 具 有 较 高 的 时 空 分 辨 率 ， 能 够 为 进 
一 步 研究 太阳 爆发 活动 提供 高 质量 的 观测 图 像 数 据 。 本 工作 中 使 用 了 两 组 不 同时 间 段 、 不同 
= 波段 、 不 同 观测 位 置 的 日 冕 喷 流 实验 数据 集 ; 
一 数据 集 I : 波段 为 171A 和 193A 两 个 波段 , 视 场 范围 449 x 179 BI, 像 元 分 辨 紊 为 0.6 
rs arcsec/pixel， 视 场 中 心 坐标 (205.2 arcsec, -253.8 arcsec)。171A 波 段 的 观测 时 间 为 2014/10/04 
10:10:00(UT) 到 2014/10/04 10:41:00 (UT), FE 156 帧 图 像 ，193A 波 段 观 测 时 间 为 2014/10/04 
10:10:08(UT) 到 2014/10/04 10:41:08 (UT)， 共 156 帧 图 像 ， 时 间 分 辩 率 都 为 12s。 数 据 集 工 : 
波段 为 171A 和 193A 两 个 波段 ， 视 场 范 围 186 x 313 像 元 ， 像 元 分 辩 率 为 0.6 arcsec/pixel, 
视 场 中 心 坐 标 (-219 arcsec, -366.9 arcsec)。171A 波 段 的 观测 时 间 为 2014/03/03 19:10:00(UT) 
到 2014/03/03 20:05:00 (UT), 4 276 帧 图 像 ，193A 波 段 观 测 时 间 为 2014/03/03 19:09:56(UT) 
到 2014/03/03 20:04:56 (UT), $ 276 帧 图 像 ， 时 间 分 辩 率 都 为 12s。 
这 两 组 数据 中 喷 流 变化 尺度 较 大 、 像素 强度 增加 比较 明显 , 同时 还 存在 一 些 亮度 增强 的 
非 日 冕 喷 流 变化 ， 非 喷 流 的 日 晃 变 化 会 导致 视 场 范围 内 的 图 像 强度 发 生 复杂 变化 。 图 1(a) 
显示 了 数据 集 1 中 在 2014/10/04 10:25:00(UT) 时 刻 171A 波 段 对 数 对 比 度 增强 后 的 全 日 面 图 
像 ; (b) 显 示 了 (a) 中 和 矩形 区 域 视 场 范围 内 171A 波 段 在 该 时 刻 对 应 的 对 数 对 比 度 增强 后 的 日 冕 
图 像 ，(c) 则 显示 的 是 193A 波 段 在 2014/10/04 10:24:56(UT) 时 刻 对 应 的 对 数 对 比 度 增 强 后 的 
日 冕 图 像 ，(b) 和 (cj 中 的 日 晃 喷 流 在 像素 强度 上 分 布 不 均匀 ， 且 (b) 中 的 图 像 内容 更 复杂 ， 日 
冕 喷 流 与 图 像 背景 像素 强度 差异 大 ， 易 于 区 别 ; 在 10:27:20(UT) 时 193A 波 段 中 日 冕 喷 流 强 
度 较 弱 ， 与 图 像 背 景 强度 差异 小 ， 此 时 171A 波 段 中 的 日 晃 喷 流 像 素 强度 与 图 像 背 景 依 然 有 
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较 大 的 区 别 。 图 2 显示 的 是 数据 集 卫 所 对 应 的 日 狗 图 像 ，(b) 中 的 图 像 内 容 比 (c) 中 的 复杂 ， 
但 (c) 中 的 日 园 喷 流 在 直观 上 易于 辨别 。 


1 (a): 数据 集 I 在 10:25:00(UT) 时 刻 171A 波 段 的 全 图 2 (a): Boi TE 19:31:49(UT) 时 刻 171A 波 段 的 全 日 面 
ARAMA; (b): 为 (a) 中 红色 矩形 框 标记 的 日 冕 图 像 ; (0): KVR: (b): 为 (3) 中 红色 和 托 形 框 标记 的 日 冕 图 像 ，(o): 为 
为 (3) 中 红色 拢 形 框 在 10:24:56(UT) 时 刻 193A 波 段 标记 的 (a) 中 红色 矩形 框 在 19:31:56(UT) 时 刻 193A 波 段 标记 的 日 冕 
HR 图 像 
Fig.1 (a): Full-disk 171A image of Data I at 10:25:00 UT. (b): Fig.2 (a): Full-disk 171A image of Data II at 19:31:49 UT. (b): 
171A image in the red rectangular box in Fig.1 (a). (c): 193A 171A image in the red rectangular box in Fig.2 (a). (c): 193A 
image at 10:24:56 UT in the red rectangular box in Fig.1 (a). image at 19:31:56 UT in the red rectangular box in Fig.2 (a) 


2 快速 鲁 棒 性 主 成 分 分 析 方法 


RPCAD9 方 法 的 数学 模型 如 式 (TD) 所 示 ; 
argmin rank(L) 4 AISI, st. D=L+S (1) 
LS 


(1) 式 是 一 个 非 凸 的 优化 问题 ， 针 对 该 问题 目前 没有 有 效 的 解决 方案 。 可 以 通过 松弛 技 
AR, KO) SPER SAL, 范 数 使 用 7 范 数 代 蔡 , 使 用 核 范 数目 ,代替 rank) KSEI 
阵 工 的 秩 ， 由 此 可 以 将 (1) 式 松弛 为 一 个 凸 优化 问题 (2) 进行 求解 。 


arg min Ii. +A IST, st. D=L+8 (2) 
L,S 
D 为 观测 矩阵 , | 如 ,是 矩阵 工 的 核 范 数 , |S 是 矩阵 5 WIL 范 数 ,4(4 > 0) 为 权重 参数 ， 


用 于 调节 低 秩 和 稀疏 的 占 比 。 其 假设 为 给 定 的 m 个 数据 ， 每 一 个 数据 列 向 量化 为 矩阵 
De R”" 的 一 列 , ERE D FY LAG) WATER PR EE L AIRES, BU D=L+S. HAL, S 
ARE ARAN, L 要 满足 低 秩 的 特性 ，$ 要 满足 稀疏 特性 ， 其 数值 可 以 为 任意 大 小 ， 这 也 是 
RPCA 和 鲁 棒 性 的 体现 。 

由 于 RPCA 是 通过 代价 高 昂 的 优化 来 计算 的 ， 在 面 对 诸 如 视频 等 数据 量 比较 大 的 情况 
时 ， 其 实时 性 比较 低 ， 因 此 RPCA 的 快速 算法 的 开发 就 显得 比较 重要 。 文 山中 对 (2) 式 的 无 
罚 约束 进行 改进 ， 得 到 在 视频 背景 建 模 /前 景 检 测 下 比 IALM 快 一 个 数量 级 的 Fast RPCA Fi 
法 ， 其 数学 模型 如 式 (3) 所 示 。 


arg min ZL +S- D], 4 AISI, st 
LS 
由 于 性 能 的 原因 ， 文 中 主要 考虑 rank(L) =t 时 的 情况 ， 如 式 (4) 所 示 。 
arg min 了 | +S- DI, AISI, strank(L)=t (4) 
L,S 


对 式 (4) 使 用 交 蔡 最 小 化 方式 进行 求解 ， 具 体 求解 过 程 如 式 (5) 式 (6) 所 示 : 
Ln = argmin|L 4S, — DI), st. rank(L)=t (5) 
L 


„St (3) 
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Sa = arg min |L, +S—D |, + AISI, (6) 


式 (5) 通 过 对 D S, BATIN SE AMR KEREL , t IAEA RAN, 
算法 有 着 非常 高 的 计算 效率 。 式 (6) 通 过 对 D-L, BET RR AY ew HR He 
(shrink(D—L,.,,) yk SBE S 。 软 阔 值 函数 定义 如 下 


di (7) 

虽然 Fast RPCA 的 实时 性 比较 高 , 但 是 ! 值 和 ,4 信 的 设置 对 实验 结果 显得 尤为 
当 <10* 时 或 者 近代 次 数 > 300 时 个 上 检测 

F 

Fast RPCA 方法 是 根据 低 秩 和 稀 琉 分 解 的 思想 来 进行 运动 物体 的 检测 ， 其 基本 思路 是 假 
设 由 序列 图 像 所 组 成 的 观测 矩阵 万 可 以 被 分 解 为 低 秩 背 景 矩 阵 志 和 稀 琉 前 景 算 阵 S 。 由 于 
太阳 日 时 时 刻 都 在 发 生 复杂 随机 的 变化 , 所 以 日 园 图 像 中 的 内 容 是 复杂 多 变 的 .从 数据 集 1 、 
I[ 中 的 序列 图 像 来 看 , 日 晃 喷 流 在 时 间 序列 上 变化 尺度 大 ， 而 序列 图 像 中 日 园 变 化 尺度 比较 
小 的 部 分 在 时 间 序列 上 差异 不 大 ， 且 在 序列 图 像 有 着 比较 高 的 占 比 , 因此 变化 尺度 比较 小 的 
这 一 部 分 日 园 内 容 可 以 近似 认为 低 秩 、 相关 的 。 而 日 园 喷 流 序列 图 像 数据 的 这 一 特点 与 Fast 
RPCA 的 低 秩 、 稀 琉 分 解 思想 相似 ， 因 此 可 以 尝试 将 Fast RPCA 方法 用 于 日 园 喷 流 的 检测 ， 
其 方法 的 可 行 性 和 优势 ， 将 在 后 续 实验 的 对 比分 析 中 进行 验证 。 


shrink(x,€) = sign(x) max {0, x 
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在 使 用 Fast RPCA 算法 检测 日 哆 喷 流 运动 之 前 ， 要 先 对 Level 1.0 的 原始 图 像 数据 进行 
预 处 理 。 预 处 理 包 括 Fits 文件 的 标准 化 、 太 阳 自 传 校正 、 曝 光 时 间 校 正和 图 像 对 比 度 增强 。 
对 于 Fits 文件 的 标准 化 和 太阳 自传 校正 ， 太 阳 软 件 包 中 提供 了 相应 的 解决 方法 (aia_prep.pro 
和 drot_map.pro)。Fits 文件 标准 化 主要 是 实现 fits 文件 日 面 中 心 对 齐 ， 像 元 分 辩 率 统一 为 0.6 
arcsec 和 消除 侧 倾角 。 太 阳 自 传 校 正 ， 主 要 是 消除 由 太阳 自传 在 投影 平面 造成 的 位 移 ， 减 小 
对 日 冕 喷 流 演化 分 析 时 的 影响 。Fits 头 文 件 中 包含 的 图 像 数据 曝光 时 间 大 小 不 一 定 相 同 ， 这 
导致 图 像 序列 的 亮度 差异 大 , 把 图 像 数据 除 以 曝光 时 间 得 到 单位 时 间 内 的 曝光 量 的 图 像 数 据 ， 
然后 在 使 用 标准 差 归 一 化 方法 ， 使 得 图 像 序列 的 强度 统一 。 原 始 fits 文件 的 图 像 数 据 的 对 比 


© 度 比 较 低 , 这 里 使 用 对 数 函 数 对 图 像 进行 对 比 度 增强 , 方便 进行 图 像 内 容 的 观察 和 喷 流 的 检 
测 。 


对 图 像 进行 预 处 理 之 后 ,使 用 Fast RPCA 算法 对 序列 图 像 进行 喷 流 检测 。 使 用 Fast RPCA 
算法 分 别 对 数据 集 [ 、 开 进行 实验 ,对 检测 结果 进行 BM3DU2 去 噪 后 ， 使 用 OTSUS HE 
像 ， 并 将 二 值 化 后 连通 域 面 积 小 于 20 IR, 得 到 最 终 检 测 结果 ,如 图 3、 图 4 所 示 。 
本 文 以 两 组 不 同时 间 段 、 不 同 波段 、 不 同 观测 位 置 的 日 冕 图 像 序列 作为 检测 对 象 ， 用 
Fast RPCA 的 检测 结果 与 帧 间 差 分 法 的 结果 进行 对 比分 析 , 来 说 明 Fast RPCA 算法 的 有 效 性 
和 优势 。 从 图 3、 图 4 中 第 二 行 Ground truth image( 使 用 MATLAB 的 Image labeler 应 用 程 
序 中 的 像素 标记 工具 ， 将 有 喷 流 的 地 方 直接 标记 出 来 ， 标 记 的 地 方 值 为 1， 其 余地 方 为 0， 
将 标记 结果 作为 标签 使 用 ) 和 第 三 行 (Fast RPCA 检测 结果 ) 图 像 序列 对 比 来 看 ，Fast RPCA 能 
够 检测 出 喷 流 的 绝 大 部 分 ， 且 与 帧 间 差 分 法 (第 四 行 ， 当 前 时 间 帧 与 前 12 秒 相 减 ) 检 测 出 的 
喷 流 区 域 相 比 ， 空 洞 较 少 ， 且 检测 出 的 区 域 比 较 连 续 ， 说 明 RPCA 能 够 用 来 检测 喷 流 。 
对 于 数据 集 1, ME 3 左 图 中 的 矩形 框 标记 的 日 园 喷 流 的 检测 结果 来 看 ， 在 10:27:12 
时 刻 ，Fast RPCA 的 检测 结果 明显 优 于 差分 法 的 检测 结果 且 检 测 结果 中 离散 点 较 少 ， 检 测 出 
的 喷 流 区 域 更 加 完整 ; 在 图 3 右 图 中 的 矩形 框 标记 的 日 园 喷 流 的 检测 结果 显示 ，Fast RPCA 
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图 3 数据 集 T ， 左 ; 171A 波 段 检测 对 比 结果 ; A: 193A 波 段 检测 对 比 结果 


Fig.3 Detection and comparison results of Data | . Left: 171A images. Right: 193A images. 

此 时 能 够 检测 出 喷 流 的 微小 区 域 ， 差 分 法 则 没 能 够 检测 出 。 对 于 数据 集 工 ， 从 图 4 REI 
框 标记 的 日 冕 喷 流 的 检测 结果 来 看 ， 其 反应 的 情况 与 图 3 矩形 框 标 记 的 日 园 喷 流 的 检测 结 
果 基 本 一 致 。 从 图 4 上 图 中 和 珑 形 框 标记 的 日 铣 喷 流 的 检测 结果 来 看 , Fast RPCA 能 够 检测 出 
更 多 的 喷 流 ,其 检测 结果 中 不 会 存在 大 量 的 非 喷 流 区 域 。 在 图 4 下 图 中 的 和 矩形 框 标记 的 日 园 
喷 流 的 检测 结果 显示 ,此 时 Fast RPCA 与 差分 法 相 比 能 够 检测 出 日 冕 喷 流 ,， 且 此 时 的 日 冕 喷 
训 检 测 结果 从 直观 上 容易 辨认 ， 不 存在 误 检 的 情况 。 
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图 4 BOSE IL, DB: 171A 波 段 检测 对 比 结果 ; 底部 ，193A 波 段 检测 对 比 结果 


Fig.4 Detection and comparison results of Data II. Top: 171A images. Bottom: 193A images. 


从 图 3. A 4 Hani mi KI RRA, Fast RPCA PEMER 28 & AE A eh 
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原理 上 讲 ， 帧 间 差 分 法 在 面 对 图 像 内 容 
一 、 背 景 稳定 的 图 像 数据 时 具有 较 好 的 


m (ah E 


# 1 Fast RPCA 与 帧 间 差 分 检测 结果 综合 比较 
Tab.1 Comprehensive comparison of Fast 
RPCA and running difference detection results 


用 性 , 但 是 在 面 对 包含 复杂 随机 变化 的 数据 集 波段 方法 Average 
= = Sas ih AS ULE YA 
SAN, 却 存在 漏 检 (图 4 E FRERE PRE1 RET Fit 
4) 的 情况 。 RPCA 0.7480 0.2514 0.3677 
en ‘ , 171 
但 是 Fast RPCA 和 帧 间 差 分 法 的 检测 差分 法 0.3939 0.3292 01585 
a i » I 
“ag 相差 不 大 (如 wi 1 RPCA 0.7596 0.1998 0.2832 
Al 4 上 下 图 (a)、(b))， 差分 法 在 检测 出 des 
Bi WS eg : ` pi 差分 法 0.6259 0.1863 0.2410 
喷 流 区 域 的 同时 ， 也 会 检测 出 许多 不 是 喷 
流 的 区 域 。 为 了 对 Fast RPCA 和 帧 间 差 分 ee oe 0.7414 0.3955 0.4877 
法 的 检测 结果 进行 客观 分 析 , 同时 也 为 Fast i 差分 法 0.3075 0.4404 0.1658 
RPCA 算法 中 的 1、 的 设置 提供 客观 依据 103 RPCA 0.6924 0.4090 0.3390 
RA E HERRN AE E ARET KA 差分 法 0.4558 0.4470 02014 
评估 。 公 式 如 (8) 式 、(9) 式 所 示 : 
TP TP 
RE = ,RE = (8) 
TP+ EP TP+ FN 
PRE. RE 
IE (9) 
PRE+RE 


其 中 ， PRE 为 准确 率 (precision)，RE AA EX (recall), FN 


综合 性 能 指标 。 7P (true 


positives), FP (false positives) 和 FN (false negatives) 分 别 表示 被 正确 检测 为 日 冕 喷 流 的 前 景 


像素 点 ， 被 错误 检测 为 日 里 喷 流 的 背景 像素 点 、 被 错误 检 涡 


1 为 背景 


的 前 景 像素 点 。 当 综合 性 


能 指标 五 越 高 ， 表 明 所 使 用 方法 的 检测 效果 越 好 。 如 表 1 所 示 。 对 于 数据 集 Il, Fast 


RPCA 在 综合 性 能 指标 上 均 有 提高 , 且 日 冕 喷 流 检测 的 准确 


验 ， 根 据 石 指标 高 低 设 置 Fast RPCA 模型 中 的 !、4 值 ， 数 据 集 工 中 171A 波 段 对 应 的 ! 值 


率 ( PRE ) 明 显 提升 。 经 过 多 次 实 


为 30, A 为 0.01, Si I P 171A 波 段 对 应 的 上 值 为 7,4 为 0.01。193A 波 段 都 采用 和 171A 


波段 一 样 的 值 。 


此 外 , 在 对 日 县 喷 流 进 行 统计 研究 时 ， 日 晃 喷 流 的 速度 、 长 度 和 宽度 等 参数 经 常 作 为 统 
计 对 象 , 这 里 使 用 FastRPCA 的 日 晃 喷 流 检测 结果 的 二 值 图 像 的 最 小 面积 外 接 矩 形 估算 喷 流 


的 长 宽 比 ， 继 而 估计 喷 流 的 大 致 形态 。 最 后 根据 检测 结果 日 
十 算 吗 流 的 总 体 速度 。 


= 


HER 


根据 估算 结果 得 到 结论 : HERMIR E KARA TE 1, 


thre i aN, 


多 态 上 具有 较 明 显 的 朝向 性 。 在 


数据 集 工 的 171A 波 段 测 得 的 速度 为 174 kmy/s， 长 宽 比 为 3.19， 在 10:25:00 时 达到 最 大 长 度 
约 85 Mm， 活 动 时 长 为 10:14:00 到 10:28:00; 在 193A 波 段 测 得 的 速度 为 194 km/s， 长 宽 比 
为 3.25， 在 19:33:44 时 达到 最 大 长 度 约 97 Mm， 活 动 时 长 为 10:13:31 到 10:26:03 。 在 数据 


集 开 的 171A 波 段 测 得 的 速度 为 174 km/s， 长 宽 比 为 4.74， 在 19:33:12 时 达到 最 大 长 度 约 93 
Mm, 活动 时 长 为 19:25:24 到 19:34:48; 在 193A 波 段 测 得 的 速度 为 157 km/s, 


发 宽 比 为 3.15， 


在 19:32:20 时 达到 最 大 长 度 约 80.3 Mm， 活 动 时 长 为 19:22:44 到 19:34:44。 


4 结论 
ASCH Fast RPCA 应 用 于 日 冕 喷 流 检测 ,通过 低 秩 和 稀 疏 分 解 的 思想 对 日 冕 变化 的 动态 


背景 进行 建 模 ， 从 而 检测 出 作为 前 景 变化 的 日 冕 喷 流 。 实 验 结果 表明 Fast RPCA 在 


有 着 喷 流 


变化 尺度 较 大 、 像素 强度 增加 明显 , 同时 还 存在 亮度 增强 的 小 尺度 非 喷 流 变 化 的 日 蝎 序 列 图 


像 中 检测 日 


冕 喷 流 时 ， 其 检测 结果 与 帧 间 差 分 法 相 比 ， 


域 连续 ， 


念 测 出 的 日 蛙 喷 流 区 域 孔洞 较 少 、 区 


日 检测 的 准确 度 有 所 提高 。 差 分 法 在 原理 上 简单 、 易 于 实现 ， 但 是 在 针对 日 冕 喷 流 


检测 这 种 除了 日 冕 喷 流 的 变化 之 外 , 还 存在 其 他 小 的 、 亮 度 增 加 的 非 喷 流 变化 的 动态 场景 时 
效果 不 可 观 ， 而 Fast RPCA 却 能 够 检测 出 帧 间 差 分 法 不 能 检测 到 的 强度 较 弱 的 日 冕 喷 流 。 


Fast RPCA 目前 存在 的 问题 是 模型 中 ! 值 、 


4 值 参数 的 设置 ， 这 些 参数 还 不 能 够 实现 自 


动 设置 ,需要 通过 多 次 试验 人 工 选 择 参数 ,而 如 何 通过 客观 指标 分 析 对 算法 的 参数 进行 自动 
设置 和 算法 的 可 行 性 进行 客观 评价 是 算法 改进 的 重点 。 而且, 日 冕 喷 流 的 提取 还 需要 进一步 


结合 日 锡 喷 流 形态 和 灰 度 变 化 的 特点 , 将 喷 流 从 检查 结果 中 完整 的 提取 出 来 ;日冕 喷 流 的 观 


SEALS pcs 


测 特征 参数 (如 : 速度 、 长 度 和 宽度 等 参数 ) 还 有 需要 更 加 精准 的 测量 。 此 外 ，Fast RPCA 在 
数据 量 比较 大 的 情况 下 运算 量 大 ,耗费 时 间 长 ,而 如 何 对 全 日 面 上 的 日 园 喷 流 ,特别 是 弱 的 
行 检测 ， 将 作为 后 续 进 一 步 的 工作 内 容 。 
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Abstract: In this paper, fast robust principal component analysis (Fast RPCA) is used to 
detect coronal jet activity in coronal sequence images. The basic idea of detection is to combine 
the idea of low rank and sparse decomposition in Fast RPCA method with the characteristics of 
coronal sequence images, such as the random background component with smaller scale and 
larger proportion, and the coronal jet with larger scale and smaller proportion, so as to realize the 
separation between the random and complex dynamic background and sparse moving objects, So 
as to detect the coronal jet as the foreground change. Two sets of coronal sequence images of 
different time periods, different channel and different observation positions from the atmospheric 
imaging assembly (AIA) observation equipment on the solar dynamics observatory (SDO) 
satellite are used as the research objects. The main research contents include the preprocessing of 
coronal image sequence, coronal jet detection, and the comparison analysis of the detection results 
between Fast RPCA method and running difference method. The experimental results show that 
compared with the running difference method, the Fast RPCA method can detect the weak coronal 
jet and improve the accuracy of coronal jet detection. 
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